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Аннотация 

Введение. Математические модели и методы повсеместно используются для исследования природных объек- 
тов, заменяя более дорогие натурные эксперименты. Одними из трудностей, возникающих при моделировании 
процессов в сложных системах, являются наличие входных данных и подбор параметров модели. Применение 
методов усвоения данных наблюдений является одним из способов оснащения математических моделей вход- 
ными данными и значениями параметров. Цель настоящего исследования состоит в прогнозировании на основе 
методов математического моделирования развития сложных природных систем в условиях загрязнения вред- 
ными веществами. Для достижения цели были решены следующие задачи: выбран метод усвоения данных 
наблюдений, актуализирована математическая модель биологической кинетики, данная модель скомплексиро- 
вана с моделью гидродинамики, разработан программный комплекс. Актуальность работы заключается в при- 
менении нового подхода к реализации модели динамики фитопланктонных популяций (эвтрофикации) Азов- 
ского моря при наличии загрязняющих примесей, основанного на применении вариационных методов усвоения 
данных, полученных в ходе экспедиционных исследований. 

Материалы и методы. Распространение загрязняющих веществ моделируется на основе трехмерной матема- 
тической модели, основанной на системе уравнений конвекции — диффузии — реакции. На входе модели по- 
дается вектор движения водной среды. Составляющие вектора скорости течений в прибрежной системе рассчи- 
тываются на основе математической модели гидродинамики, базирующейся на трех уравнениях движения и 
уравнении неразрывности. Разработанный на основе описанных моделей программный комплекс получает на 
входе натурные данные, собранные в ходе экспедиционных исследований, и позволяет уточнять модель загряз- 
нения водной среды и биоты благодаря применению вариационных методов усвоения данных. 

Результаты исследования. Построен краткосрочный прогноз распространения загрязняющих веществ на вы- 
ходе из Таганрогского залива. Проведенный вычислительный эксперимент отражает динамику распростране- 
ния загрязняющих веществ от источников заражения на временном интервале от 3 до 12 дней. 

Обсуждение и заключение. Рассмотренные в данном исследовании вариационные методы усвоения данных 
наблюдений позволяют уточнять и дополнять математические модели динамики фитопланктонных популяций 
и распространения загрязняющих веществ. Программное обеспечение, основанное на описанных в данной ра- 
боте математических моделях, дает возможность строить кратко- и среднесрочные прогнозы распространения 
вредных примесей, оценивать их влияние на развитие основных видов фитопланктонных популяций в Азов- 
ском море и определять стратегии управления устойчивым развитием. 


Ключевые слова: модель эвтрофикации, модель гидродинамики, вариационные методы, опасные явления, 
усвоение данных наблюдений 
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Введение. Математические модели и методы уже несколько десятилетий успешно применяются для проведе- 
НИЯ исследований в различных сферах науки и инженерии. Математическое моделирование предоставляет быст- 
рый, удобный и относительно недорогой инструментарий для изучения и прогнозирования процессов, протекаю- 
щих в сложных, по сравнению с экспедициями и натурными экспериментами, природных системах, и применяет- 
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ся для решения многих научных и практических задач. Так, например, могут быть решены задачи прогностиче- 
ского моделирования заиления судоходных путей, что важно для безопасного судоходства, а также предсказания 
последствий чрезвычайных ситуаций и катастроф техногенного характера. Примером может служить сильный 
шторм 11 ноября 2007 года в азово-черноморском бассейне, в результате которого более 20 судов потерпели кру- 
шение, а район Керченского пролива стал местом экологической катастрофы. В воду попало несколько тонн мазу- 
та и серы, в результате чего произошло заражение береговой линии и седиментного слоя соединениями нефтепро- 
дуктов, последствия его наблюдались еще несколько лет. Другой пример неблагоприятных процессов — транс- 
порт донных материалов из устья реки Дон в Таганрогский залив, что влияет на ареалы гидробионтов, способ- 
ствует интенсивной эвтрофикации и приводит к размножению комара-звонца (СЫтопопиае Ме\лтап). 

Согласно Постановлению Правительства РФ № 2451 от 31.12.2020 г.!, всего лишь четыре часа с момента 
обнаружения или с момента поступления информации о разливе в водоеме есть у исследователей, ответствен- 
ных лиц, принимающих решения, а также представителей водоохранных служб для расчета изменения концен- 
трации нефти и нефтепродуктов и их локализации. Согласно приведенному и другим нормативным докумен- 
там, принятым Правительством РФ, ответственные лица и структуры должны принять решение и произвести 
действия по устранению опасной экологической ситуации природного и техногенного характера в течение не- 
скольких часов-суток. Исходя из этого время построения прогнозов и сценариев развития чрезвычайной ситуа- 
ции (ЧС) ограничено. Это требование определяет актуальность разработки комплекса математических моделей 
гидродинамики и гидробиологии, которые учитывали бы особенности прибрежных систем (воздействие ветров 
на структуру течений, силу Кориолиса, сложную геометрию расчетной области, турбулентный обмен, испаре- 
ние, сгонно-нагонные явления, стоки рек и др.) и позволяли получать прогнозы за минимальное время. 

При построении прогнозов развития природных систем и попытке описать реальное физическое явление 
средствами математического моделирования зачастую только лишь построения функции состояния моделей 
процессов недостаточно. Повысить точность решения позволяет использование для математических моделей 
сопряженных задач, а также алгоритмов, базирующихся на вариационных принципах. Эти принципы позволя- 
ют установить связь модели с натурными данными [1]. Такая методология является эффективной при решении 
прикладных, в том числе и вычислительно трудоемких задач, например вариационных и оптимизационных за- 
дач математической и ядерной физики [2]. Г.И. Марчуком и его последователями применялись сопряженные 
уравнения [3], что позволило повысить эффективность решения задач аэро- и гидрофизики для атмосферы и 
глубоководных водоемов, а также была усовершенствована теория построения сопряженных операторов для 
линейных и нелинейных моделей [4]. 

Вариационный подход к решению объединенных прямых и сопряженных задач с использованием методов 
усвоения позволил улучшить связь между математическими моделями и натурными данными. Методы усвое- 
ния данных, которые развиваются с 1960-х годов, основаны на построении обратных и оптимизационных задач 
с использованием двух подходов: классического вариационного принципа Лагранжа с применением сопряжен- 
ных задач [4] и оптимизационных методов типа взвешенных наименьших квадратов [5]. 

Усвоение данных наблюдений является инструментарием, который позволяет существенно повысить точ- 
ность прогностического моделирования природных процессов, он давно и успешно применяется ученым сооб- 
ществом [6]. Здесь актуальной задачей становится разработка новых методов, которые позволили бы суще- 
ственно сократить время расчетов. 

Материалы и методы. При построении моделей прогнозирования природных явлений и процессов одними 
из основных проблем являются вопросы соответствия решения, полученного с использованием математической 
модели, реальному процессу, протекающему в природной системе, и уменьшения процента неопределенностей. 

При построении математических моделей гидродинамических и гидробиологических процессов требуется 
информация о начальных условиях и параметрах модели, которая может быть получена с помощью данных 
наблюдений. Таким образом, при построении прогностических сценариев протекания ЧС природного или тех- 
ногенного характера очень важно оценить адекватность самой математической модели. Следующий этап моде- 
лирования включает в себя проверку корректности, устойчивости поставленной задачи. Исследование матема- 
тической модели на непрерывном уровне предполагает исследование влияния входных данных на решение мо- 
дельной задачи. Возмущение правых частей используемого уравнения или системы уравнений в частных про- 
изводных в рассматриваемой задаче Коши позволяет при известных операторах исследовать свойства постро- 
енной математической модели. Исследование стационарных и особых точек непрерывной функции или 
нескольких функций — решений поставленной задачи, например концентрации одного или нескольких загряз- 
няющих водную среду веществ, позволяет разработать сценарии, от пессимистичного до оптимистичного ре- 
жимов, с целью разработки мер эффективного управления сложной водной экосистемой. 


1 Об утверждении Правил организации мероприятий по предупреждению и ликвидации разливов нефти и нефтепродуктов на территории 
Российской Федерации, за исключением внутренних морских вод Российской Федерации и территориального моря Российской Федерации. 
Постановление Правительства РФ № 2451 от 31.12.2020 г. ОВГ: В рз://40с$.сп.га/4осиаеп 573319208 (дата обращения: 21.05.2024). 
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Впервые для анализа натурных данных был применен метод полиномиальной интерполяции данных посто- 
янно пополняемой базы экологических замеров в двумерном случае. Для наблюдений были определены обла- 
сти влияния. При реализации алгоритма для расчета в текущий момент БасКэгоипа — данные ранее полученно- 
го прогноза — использовались как входная информация. 

Методология решения задач усвоения данных вышла на новый уровень с появлением метода ОГ (Орйта! ш- 
{егро!айоп), или метода статистической интерполяции. 

Следующий этап развития рассматриваемых методов связан с развитием и реализацией вариационных ме- 
тодов, включая теорию оптимального управления. Эти методы базируются на минимизации функционала, по- 
строенного специальным образом, с помощью которого устанавливается связь между решениями и наблюдени- 
ями (натурными замерами, экспедиционными данными, базами геоинформационных систем ГИС). 

Данная теория и методы широко используются при реализации задач метеорологии [5] и динамической оке- 
анографии [7]. В процессе минимизации построенного функционала необходимо вычислять его градиент, для 
чего успешно применяются сопряженные уравнения, что описано в работах [8, 9]. 

Модель динамики фитопланктонных популяций. Математическая модель динамики фитопланктонных 
популяций описывает процесс активного роста микроводорослей при наличии достаточного количества био- 
генных элементов. Если развитие фитопланктонных популяций становится слишком интенсивным, говорят о 
процессе эвтрофикации. Причины эвтрофикации могут иметь как природный (климатические изменения), так и 
антропогенный характер (поступление в водоем значительного количества биогенных веществ со стоками рек). 
В системе С; — значения концентрации 1-ой субстанции [10, 11]: 


д + С 
9Сь , 0(иСь) (УС, ((®+ус»)Сь) д (= гс), 8[, 2), д |. гс.) " 
[911 бх ду 02 бх бх ду ду 02 02 


где и = {и, у, и} — вектор скорости среды (водного потока); исх — гравитационное осаждение А-ой компонен- 
ты, если она находится во взвешенном состоянии; 1, у, — горизонтальная и вертикальная составляющие коэф- 


фициента турбулентного обмена для К-ой компоненты; \ух — химико-биологический источник (сток) или член, 
описывающий агрегирование (слипание-разлипание), если соответствующая компонента является взвесью, ин- 
декс А указывает на вид субстанции, К = 1,15: 

1 — сероводород (Н25); 

2 — элементная сера 5; 

3 — тиосульфаты (и сульфиты); 

4 — сульфаты (504); 

5 — общий органический азот (№); 

6 — аммоний (№Н4 — аммонийный азот); 

7 — нитраты (№Оз); 

8 — нитриты (№О>); 

9 — фитопланктон; 

10 — зоопланктон; 

11] — силикаты (51Юз — метасиликат; 5104 — ортосиликат); 

12 — растворенный кислород (05); 

13 — железо (Ёе?*); 

14 — фосфаты (РО4); 

15 — кремнекислота (Н>5 04, Н›5Юз — ортокремневая и метакремневая кислоты соответственно); 

16 — микропластик. 

Система (1) содержит уравнения, которые можно отнести к типу конвекции — диффузии — реакции. В качестве 
расчетной области рассмотрим замкнутый бассейн С. Невозмущенная поверхность водоема Ху ограничивает сверху С, 
Уь = Ух, у) — донная поверхность снизу. о — цилиндрическая поверхность, ограничивает С’ сбоку. Введем обозначе- 
ние: Х = ХЦ с Ц УХ» — кусочно-гладкая граница области С’, временной интервал 0 <{< То. Полагаем, что п и ии — век- 
тор внешней нормали и нормальная составляющая и к поверхности У. 


Считаем, что начальные условия для системы (1) выглядят следующим образом: Сино = Сю(х, у, 2), К = 1,15. 


Скомплексируем (1) со следующими комбинированными граничными условиями: на о: если < 0, то 


Сь= 0; если ив > 0, то 0; на У: Св (Сь); на дне Хы: РЕ в, ‚К =1,16, определим =х как коэф- 
и = 2 


фициент поглощения К-ой примеси донными отложениями. 
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При безветрии, особенно в летний период, в придонных слоях мелководных водоемов, таких как, например, 
Азовское море, Таганрогский залив, Геленджикская бухта, могут возникнуть практически анаэробные условия. 
Восстановление водонасыщенного поверхностного ила влечет за собой высвобождение в раствор железа, фос- 
фатов, сульфатов, марганца, аммония и силикатов, а также органических соединений. Скомплексированные 
модели вида (1) (разработанная авторским коллективом и модель гидродинамики [12]) используются для изу- 
чения механизмов окисления и восстановления марганца, ассимиляции №МНа, нитрификации, нитратредукции 
(денитрификации), аммонификации, окисления Н>5, сульфатредукции и др. Эксперименты с моделью (1) дают 
возможность изучать биогенный и кислородный режимы прибрежной системы, анализировать механизм фор- 
мирования вследствие антропогенной эвтрофикации заморных явлений рыб и других гидробионтов. 

Вариационный подход для пространственно-трёхмерной математической модели эвтрофикации вод. 
Запишем математическую модель (1) для расчетной области (Азовское море) в виде операторного уравнения: 


Исх) = рен х)-у-г=0, ©) 


определим С как вектор-функцию состояния изучаемой водной экосистемы С= [с р (х,1),К = 116} ь 
С=С(х,!) Е О(Ц,), Ц, =Сх(0,Ть), (х,!)еЦ,; /($,У) — дифференциальный нелинейный пространствен- 
ный оператор; р — диагональная матрица; у = [у й (х,г),К = 116} — вектор, компонентами которого являются 


нкции источников; Г=уИ, (Х,1 ‚К = 1,16 —Щ вектор, компоненты которого содержат нкции неопределенно- 
> К 


стей и ошибок математической модели (1) с начальными и граничными условиями. Зависимости (модели 
наблюдений), коэффициенты и параметры ис», щ, и, у, и’, уг входные данные начальных и краевых условий ДлЯ 


модели (2), К = 1,16; внутренние параметры операторов включаютсяв УеК(Ц,). 
Пусть 1 = 0, тогда начальные условия для (2) будут иметь вид: 
С = С + У=У2 +6, (3) 
где С’и У* — априорные оценки вектор-функции состояния и вектора параметров соответственно; функции 


неопределенностей обозначены с помощью &ё и 6. 
Рассмотрим интегральное тождество: 


Ксх,С*) = [(ссу). сада =0, (4) 
Ц, 
здесь С” — это функции, сопряженные к С(С* Е О"(Ц))). (4) представляет собой вариационную постановку мо- 
дельной задачи (2), (3), или функционал энергетического типа. Перепишем (4) в следующем виде: 
16 
и(с.х,с")= У (лс, - [+ +. Саба = 0. (5) 
и Ц, 

Операторы турбулентного обмена и переноса входят в слагаемые (ЛС, С"). 

Модель гидродинамики [12] будем считать моделью процесса. В моделях гидробиологии параметризацию 
параметров и в результате полученные функциональные зависимости, например, для описания продукционно- 
деструкционных процессов или роста фито- или зоопланктона будем считать подмоделями, или моделями 
наблюдений. Определим зависимость между замерами и функциями состояния: 


(1 =[\(С)], +1(х,/), (6) 


где [\(С)]„ — вектор подмоделей (моделей наблюдений); Ц(х, 0) — вектор ошибок и неопределенностей; 
ф„ — величины, за которыми осуществляем наблюдение. 


т 


Определим ф»„ на Цу' ЕЦ,. В (6) операция переноса информации с Ц, на Цу’ обозначена квадратными скобками. 


Расширим систему моделирования данными натурных замеров (считаем их близкими к точным), при этом 
функционал «качества» будет иметь следующий вид: 


Ф.(©)= (© -[% (©), ми (= -[%(©)],)| = ("см ©) 
ЦР 
определим хо как весовую функцию для определения конфигурации носителя наблюдений Ц” в Ц; и интегралы 


по области Ць представляющие собой меру для (7) в виде С! = Мхо(х, 0, где М — весовая матрица. 
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Рассмотрим функционалы, представляющие собой обобщенные характеристики поведения гидробиогеоценоза: 


Ф, (С) = [5 (Суд даба = (Ежи), ие 0" (Ц,), =. 
Ц, 
ЕС) — ограниченные и дифференцируемые относительно С Е О(Ц;) функции, которые будем оценивать. 
Определим функционал для минимизации неопределенностей: 


$ (С) =Ф* (С) м -в),. +( “Сок }.+ + (< - с) (сс + 


Ци 
(8) 
+(у°-у | С. (У -У) ). /2+1' (С,У,С* к >1. 
(це) 
Считаем, что С; — весовые матрицы, {= 1,4. Рассмотрим систему: 
=рл,С+." (С,У)-у-г=0; (9) 
Ай 
29 (р, уст" (СУ) +а, =0; 
9С 
* В д т й я : 
с; (х)_‚ =0; 4, = а ® (С)+0,5(п с) 


С° = С +СУ'С; (0),1= 0; г(х,Е) = СС} (х,1); 
д : 
У=У, +СИГ,: Г, =” (С.У,С;); 


‚К>1. 


а=0 


Й д 7! 
А(С,У)С =.“ (С+аС’У)} 


АЖС, У), Л; определим как операторов сопряженной задачи и производных или их дискретных аппроксимаций 
по времени; Г, — функции чувствительности моделей к изменению параметров; С’=56С; а — заданное число. 

Рассмотрим алгоритмы усвоения данных последовательных наблюдений, поступающих из различных 
наблюдательных средств в систему моделирования в режиме реального времени. Для этого воспользуемся ме- 
тодами расщепления и декомпозиции: 


№, -—1 М,-Г р 
= де [ОЕ УФ (10) 
п=1 п= [= 

здесь Фи — = часть функционала (9) для [-ьё] на [-м этапе — расщепления, 


п=1, М, ‚ р определим как общее число этапов расщепления. Дискретизацию проведем на основе аддитивно- 
усредненных схем расщепления. Будем использовать алгоритм для нахождения решения в области Ц” с регу- 
лярной равномерной временной сеткой ©" = [т =0, № ,} . Для исследований будем использовать структуру с 


фазовыми пространствами: 
п *п й 
| к СР. ЕЫ ем г =9"щ Е) 


Метод решения поставленной задачи. Для решения (9) считаем, что п =1,№,. Опишем пошагово алго- 


ритм используемого метода. 
1. Перейти на подсеточную структуру декомпозиции, при этом =, 


[Сб 0*(Ц!}}, [ег «он #}}, ст =Ст. 


[= 
2. В подсеточной структуре получить решения прямых И сопряженных задач: 
АРС” тг =0,1=1р, р>1. 
п-1 
ик" _ | ОФи(С) 
ЛС: [| +0` с (е„- [\(<)],) , 


1 


СИН =0, 7 -(@) 62. рт. 
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п-1 


В функции у’ входят С’. г’ учитывают все неопределенности на шаге [1.1 {н|. 


3. Возврат в основную структуру О” (ц’) при 2 = &: 


©; = 0! (Ц! | > {С* 0" (ц!)}, С" тег. 
1=1 1=1 


Считаем, что последний этап расщепления можно реализовать по формуле: 
АтС" у" г" =0, а 


‘рп 


здесь ЛС" — оператор аппроксимации часть модели на р-м этапе; у, — функции источников; г, — функция 


неопределенностей модели (2) и вносимых в дискретную модель расщеплений на шаге. Считаем, что матрицы 
весов Си и Со, известны. 
Рассмотрим следующую задачу: 


рп 


АН, С 9" С} (12) 
м = (Сы а», С”. (13) 


Для решения будем использовать следующий метод, алгоритм которого опишем пошагово: 

1. Рассчитать по (12) С*". 

2. Найти г„ , пользуясь (13). 

3. Рассчитать С", используя формулу (11). 

Возникающую в процессе реализации алгоритма метода СЛАУ будем решать методом прогонки. 
Рассмотрим преобразованный алгоритм: 


АС" = аыСь (Ф"-С"), (14) 


рп 
А С" у" (Сы / а} С" =0. (15) 


Систему решим методом прогонки. Для рассмотренных модификаций устойчивость схем расщепления обеспе- 
чивает устойчивость используемых схем усвоения данных. Если в (8) С* = соп5, то оно дает балансовое соотноше- 
ние первого порядка. В случае, если С* = С, то получаем уравнение баланса энергии анализируемой системы. 

Выпишем функционал качества для новой последующей модификации: 


Фи (С) =0,5] а (пт, )+а (ик, ) |, (16) 
г, = Л мС„ —\,, 1, =Ф) =: +92 = 1 01, &2 > 0. (17) 
Найдем минимум функционала относительно функции С”: 
* п п 0 Ш п п 
АыСь (АтС" №. С», (С"-")=0. (18) 
2п 


Получим однозначно разрешимую систему с пятидиагональной матрицей. Решим полученную СЛАУ мето- 
дом прогонки. 

Новый класс методов усвоения в реальном времени включает в себя схему аддитивного последовательного 
усвоения. Построенные модификации ввиду большого объема вычислительной работы ориентированы на су- 
первычислительные системы, включая кластерные системы и графические ускорители. 

Замеры ф„ применяются в виде карт и цифровых изображений. Такое представление дает значительную 
плотность данных в области Ц ь замеры представляют собой информационные поля. Планирование наблюдений 
базируется на значениях функции неопределенностей. Там, где они большие, планируется проводить дополни- 
тельные натурные замеры или наблюдения. 

Результаты исследования. Для решения задачи моделирования эвтрофикации вод Азовского моря (1) раз- 
работан комплекс параллельных программ, включающий в себя: 

— модуль гидродинамических процессов, рассчитывающий поле течений водного потока на основе матема- 
тической модели для мелководного водоема [12]; 

— модуль распространения загрязнений в водной среде и изменения концентрации основных гидробион- 
тов (1), позволяющий оценить влияние загрязняющих веществ на биологическую продуктивность акватории; 

— карту глубин Азовского моря для построения расчетных сеток для численной реализации разработанных 
алгоритмов; 

— базу экспедиционных данных, позволяющую уточнять модель загрязнения водной среды и распростране- 
ния биоты благодаря применению описанных выше методов усвоения данных. 
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На основе разработанного ПО проведен численный эксперимент. Вектор скорости водного потока рассчи- 
тывается с помощью модели гидродинамики при восточном ветре, скорость которого составляет 5 м/с, и пода- 
ется на вход для расчета движения загрязнений, содержащих микропластиковые частицы, на основе уравнений 
конвекции-диффузии. На рис. 1 а представлены результаты численного эксперимента по расчету полей течений 
водного потока при указанных метеорологических условиях, на котором можно наблюдать вихревые структуры 
течений в районе кос, в северо-восточной части моря, в Таганрогском заливе. Градиент цвета на рис. 1 а отра- 
жает распределение скоростей водного потока, максимальное значение — 4,822 м/с. На рисунках 1 6, в, г пред- 
ставлены результаты расчета для модельной ситуации распространения заражения водной среды микропласти- 
ковыми частицами, попадающими в Азовское море со стоками рек Дон и Кубань, а также когда источник за- 
грязнения опасными веществами, в том числе микропластиком, находится на выходе из Таганрогского залива. 
Представлен краткосрочный прогноз распространения загрязняющих веществ через 3, 6, 12 суток с момента 
заражения. Начальная концентрация загрязняющего вещества составила 5 мг/л, через 3 суток максимальная 
концентрация на выходе из Таганрогского залива составила 1,363 мг/л, через 6 суток — 0,83 мг/л, че- 


рез 12 суток — 0,336 мг/л. 
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Рис. 1. Картина течений в Азовском море и распространение загрязняющих веществ, временной интервал: 
а — начальная концентрация; 6 — 3 суток; в — 6 суток; г — 12 суток 
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Обсуждение и заключение. Вычислительный эксперимент демонстрирует, что несмотря на восточный ве- 
тер, вихревые структуры течений захватили загрязняющие вещества и переместили их в Таганрогский залив. 
Устойчивые вихри потенциально могут захватывать и удерживать частицы микропластика, попадающие в море 
со стоками рек, а также способствовать накоплению загрязнений в придонном слое в результате биообрастания 
микропастиковых частиц и их затопления. 

Как отмечалось выше, при построении математических моделей для прогнозирования природных явлений и 
процессов одной из основных проблем является проверка их адекватности путем анализа полученных на их 
основе результатов на соответствие поведению изучаемой природной системы. При построении математиче- 
ских моделей гидродинамических и гидробиологических процессов требуется информация о начальных усло- 
виях и параметрах (входных данных), которая может быть получена с помощью наблюдений. Таким образом, 
при построении прогностических сценариев необходимо не только оценивать качество построенной математи- 
ческой модели, но и усваивать данные наблюдений, исследовать чувствительность построенных моделей к из- 
менениям входных данных. 

В работе представлен подход к реализации модели динамики фитопланктонных популяций (эвтрофикации) 
Азовского моря с применением вариационных методов усвоения данных, полученных в ходе экспедиционных 
исследований. Разработанный программный комплекс использует материалы экспедиционных работ, постоян- 
но пополняемые базы экологических данных, ГИС и позволяет уточнять модель загрязнения водной среды и 
распространения гидробионтов благодаря применению вариационных методов усвоения данных. Разработан- 
ный программный комплекс позволяет прогнозировать распространение загрязняющих веществ в прибрежной 
системе, некоторые из них, например биогенные вещества, являются питательной средой и способствуют раз- 
витию опасных микроводорослей. Данный прогноз позволяет вырабатывать стратегии по управлению устойчи- 
вым развитием природной системы. 
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